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Abstract 

We study the electroweak phase transition in a supersymmetric 
version of the Standard Model, in which a gauge singlet superfield is 
added to the Higgs sector. We show that the order of the transition is 
determined by the trilinear soft supersymmetry breaking terms rather 
than by the O(m^T) term in the 1-loop, T 7^ corrections. This 
fact removes the Standard Model upper bound on the Higgs mass, 
niH < 55GeV , coming from the requirement that baryon asymmetry 
is not washed out by anomalous electroweak processes. We perform a 
numerical analysis of parameter space including in the effective poten- 
tial top-stop contribution to 1-loop radiative corrections. We find that 
this model is compatible with the preservation of baryon asymmetry 
for masses of the lightest scalar up to about 170 GeV. 



1. Introduction. 

Several years ago Sakharov [1] realized that the observed baryon asymme- 
try of the Universe (BAU) could have been generated at some early stage of 
cosmic evolution at which three conditions were fulfilled: B-violating interac- 
tions, C and CP violation, and a departure from thermodynamic equilibrium. 
All these conditions are fulfilled in grand unification theories (GUT), in which 
the baryon asymmetry is generated in the out-of-equilibrium, B-violating, 
decay of some superheavy boson [2]. However this scenario presents severe 
complications. First of all, as noted by Bocharcv, Kuzmin and Shaposh- 
nikov[3], any nonzero fermion number (B+L) created at the GUT epoch is 
almost completely erased [4] by the anomalous electroweak B-|-L-violating 
processes [5], which are in equilibrium down to a temperature T ~ 10^ GeV. 
So, if no B-L is created at the GUT phase transition, then no baryon excess 
survives down to the Fermi scale. Moreover, if B-L-violating interactions 
{6L = 2 Majorana masses for neutrinos, R-parity violating interactions in 
SUSY models [5]...) are present and in equilibrium at high (T ^ 100 GeV) 
tcuipcraturc, then also an eventual B-L component of the GUT-generated 
BAU vanishes before the onset of the electroweak era. 

The above considerations led many people to investigate the possibility 
that the BAU was generated at the electroweak phase transition, ([6]-[13]). 
Here the difficulty is threefold; first, CP violation in the Standard Model is 
too small [7], second, it is not obvious that the transition is of first order 
strongly enough to make the baryon number production effective and, third, 
we must avoid the wiping out of the baryon asymmetry by the anomalous 
processes. The efficiency of the latter is suppressed by the exponential of 
(minus) the 'sphaleron' mass [6] over the temperature 

E'P^{T) _ 47rv{T) B 

T ~ gw T' ^> 

where -B(-j-) is a slowly varying function of the couphngs 5(0) = 1.56, 
B{oo) = 2.72, and v{T) the value of the vacuum expectation value (VEV) of 
the Higgs field at the temperature T. The requirement that (1) is sufficiently 
large at the phase transition gives a lower bound [6] on the value of the 
vacuum expectation value (VEV) at the critical temperature (which we define 
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here as the temperature at which the effective potential becomes ffat at the 
origin, V"{0) = 0) 



In the Standard Model v{Tc)/Tc is given by the ratio of the coefficient 
of the term cubic in the Higgs field, 6, and the quartic self-coupling, A. So it 
is easy to see that (2) imposes an upper bound on the Higgs mass 



The crucial 6 coefficient is zero in the Standard Model at the tree level; 
it is generated by the one- loop finite temperature corrections [14,15] as a 
subleading term in the M/T expansion [15], and its evaluation requires a 
careful resummation of the leading infrared (IR) divergencies [14,16]. Despite 
all the uncertainties related to this procedure, it seems by now established 
that the Standard Model bound (3) has been ruled out by LEP results, 
Mh > 57 GeV [17], or it will be in the near future. 

Since 6 receives contribution by the bosonic degrees of freedom with 
masses lighter than the temperatures in consideration 0, attempts have been 
made to relax the bound on Mh by considering the Standard Model with an 
extended Higgs sector [6,8,9,18], or minimal SUSY Standard Model (MSSM), 
in which two higgs doublets are present [10]. In all these cases the cubic term 
is again a one-loop effect and its evaluation meets the same IR problems that 
we find in the Standard Model. Moreover, as discussed in ref. [13], there is 
a large suppression (0(a^)) of the baryon asymmetry in models contain- 
ing only scalar doublets, without extra CP violation sources than the Higgs 
potential and the Yukawa couplings. 

In this paper we study the electroweak phase transition in the simplest 
possible extension of the MSSM with soft supersymmetry breaking, in which 
a gauge singlet superfield is added to the superpotential [19,20]. From the 
point of view of the phase transition the main feature of this model is the oc- 
currence of trilinear soft SUSY breaking terms in the Higgs sector of the tree 
level effective potential. These terms behave as effective cubic terms in the 
radial direction from the origin of the field space to the vacuum, so that the 

^in fact it may be shown numerically that the M/T expansion is a good approximation 
up to ^ ~ 2.2 [9] 




(3) 



3 



potential has a barrier along this direction already at the tree-level. More- 
over, including the leading (no IR problem) finite temperature corrections, at 
one-loop, we get also positive linear terms which enhance the barrier between 
the false and the true vacuum. As a result, this model leads very naturally 
to a strongly first order phase transition and, as we will see, the bound (2) 
may be fulfilled in a large portion of the parameter space corresponding to 
values of the mass of the fighter scalar up to about 170 GeV . 



2. The model. 



We assume [20] that the squark and slepton fields have vanishing vacuum 
expectation values so that we can restrict our attention to the Higgs sector 
of the superpotential involving the superfields Hi, H2 and N only 

Wmggs = XHiH2N-^kN\ 

The tree level scalar potential is given by 

V = Vp + Vn + Vsofu (4) 
Vf = \X\^\\N\^{\Hi\^ + \H2\^) + \HiH2fl+k^\N\^ - {Xk*HiH2N^* + h.cy, 



Vd = ^^^mi' - \Hi\'r + ^\HlH2\'; 
Vsoft = ml^\Hi\' + mljH2\' + m%\N\^ 

-{\AxH1H2N + h.c.) - {\kAkN^ + /i.e.), 

where//i^(^i),i/2^(^^o )• 

Note in particular the presence of trilinear terms in Vsofti which are absent 
in the Higgs sector of the MSSM. 

Redefining the global phases of Hi and we can always take XAx and 
kAk real and positive, while by an SU(2)XU(1) global rotation we put 

v+ =< H+ >= 0, V2 =< H° >e R+. 
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We will not discuss the CP violation aspect of the BAU generation prob- 
lem, so we assume CP conservation. Choosing A/c > we avoid explicit CP 
violation and find three degenerate minima (note that the problem has a 
discrete Z3 symmetry[21]) corresponding to different values for 0i and (p^, 
which are respectively the phases of Vi =< i^i > and x =< N > 

' 01 = , = 0; 

< 01 = Itt , 0x = Itt; (5) 

Due to our assumption of CP conservation, the true vacuum of the theory 
is the one with both phases equal to zero. 
We will work in the unitary gauge 

_ / (i?ei/o + -^sin/3AO) \ / C0S/3C+ 

^'~[ sin/3 C+* J ' V (^e//o + ^cos/3AO) 

where, as usual, tgP = V2/V1, and 

= V2{sm(3 ImH^ + cos (3 ImH^); 
C+ = cos/3i/+ + sin/3i/-*. 

The combinations orthogonal to and arc the would-be Goldstone 
bosons giving masses to the and The physical degrees of freedom 
of the model are given by three scalars, two pseudoscalars and one charged 
scalar. 

Before proceeding in the calculation of the finite temperature corrections, 
we wish to comment on some aspects of the tree level potential which will be 
used in what follows. First, it contains seven free parameters: A, k, Ax, Ak, 
rriff^, mjj^ and m%. Imposing the stationarity conditions in [H^^H^^N) — 
(vi, V'i, x) with the constraint 

vl + vl = v^ = (246.2 GeVf, 

we express the soft masses in terms of the six parameters A, Ax, ^fe, tg(5 
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and X 



m 



Hi 



X{Ax + kx)xtg(3 - \^{x^ + v'^ sin^ (3) 
£^^;2cos2/3 



m 



H2 



+ 



\{Ax + kx) X cotg(5 - \^{x^ + v'^ cos^ (5) 
4 v"^ cos 2(5 



(6) 



= XAx^sm2p + kAkX - X^v^ 
+ Xkv'^sm2p 

Note that (6) does not guarantee that {vi,V2,x) is the global minimum of 
the effective potential (4); for each choice of the parameters we must verify 
that this is indeed the case. 

In addition, we have to be sure that v~ vanishes in the vacuum; looking 
at the potential along the charged direction, 

,2 



ch 



\H- 



m 



o'2 



l^ + X'\N\' + ^-{\H,\' + \H,\') + 

H2\ 



+ 



9^ + 9'^ 



H- 



we see that a sufficient condition in order to have v = as a global minimum 
is 



m 



li+X'x' + ^v' + ^v'cos2p>0 



(7) 



Now we introduce the finite temperature corrections at the one-loop level. 
As usual [15] they can be split into T — and T 7^ contributions. 
The former are given by ([22]) 



Str I M\(f)) 



In 



(8) 



where M^(0), with cj) — (Hi, H2,N), is the field dependent mass matrix and 
Q the renormalization point. The dependence in (8) is compensated by 
that of the renormalized parameters, so that the full effective potential is 
independent of . 

We have approximated expression (8) by considering only the top and stop 
contributions: in addition, we take the two stops' masses to be degenerate 

ml = him', (9) 



mf mn 

tL,R V 



htm 
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where itiq 



is the squarks soft mass (we assume mfj = rriQ) and ht is the 
Yukawa couphng. We fix the renormahzation scale Q at a value Q such that 
the tree level relations (6), expressed in terms of the renormalized parameters, 
do not change, i.e. 



dHo 



H2 



_2 =0. 

Q 



(10) 



V2 



Since we have neglected the ti — mixing terms in the stop mass matrix, 
the correction (8) does not affect the pseudoscalar mass matrix [22], while it 
modifies the scalar spectrum. 

The T part of one-loop corrections [15] is given by the expression 



'>,T) 



+ 



E 

F 



E 

B 



27r2 



27r2 



dx X In 



dx X In 



1-e 



(11) 



;i2) 



where ?ti|,^^-|(0) is the tree-level mass of a fermion (boson) in presence of 
the background fields 0(= Hi, H2, N), gF(B) is the corresponding number of 
degrees of freedom and the sum runs over all fermions (F) and bosons (B) 
of the theory. For values of the fields such that m(0)/T < 1 we can expand 
(11, 12) as ([15]) 



T^o {(l>,T) ^ Y^fQfT^ 
+ Y^bQeT^ 



720 ~ 48r2 



-— + 

90 ~ 



24T2 



+ 

+ 



rp4. t I V rp 



T3 



(13) 



where we have neglected the subleading m\/T^ term in the expansion for 
the bosons. As we said, this term is essential in the Standard Model as in the 
MSSM in order to obtain a first order phase transition, while in the present 
model the barrier between the true and the false vacua is given mainly by 
the tree-level trilinear terms. We have checked numerically that the cubic 
term coming from one- loop corrections changes the values of Tc and V{Tc) 
only by a few percent. 
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The degrees of freedom corresponding to masses m > T are Boltzmann- 
suppressed and in this hmit the correction (15) reduces to ([9]) 



Since we are considering temperatures of the order of Mw the exponential 
factor in the previous expression allows us, with good approximation, to 
neglect SUSY particles with masses above the SUSY threshold which we 
assume to be 



Lighter particles will contribute to (13) even if it may happen that M(0)/T ~ 
1 for some of them, at large values of the field. However (see the footnote 
n. 1) this does not change matters essentially with respect to the situation 
in which we use (11,12). Moreover the use of (13) allows us to perform an 
analytic study of the critical temperature and of the phase transition which 
can help us to understand the main properties of the model, otherwise hidden 
by a numerical evaluation of the integrals in (11,12). 

In the following we will derive the mass matrix for the relevant degrees 
of freedom of the theory on the background of all the scalar fields. We will 
then discuss bounds on the parameter space coming from the requirement of 
a correct symmetry breaking pattern (i.e. no VEV's for A^, ImN and for 
C+), and on the LEP limits on chargino masses. Finally we will scan the 
allowed parameter space, find the values of the critical temperatures, and the 
minimum of the potential at those temperatures. 

3. The mass matrix. 



Let us consider the matrix of the second derivatives of the tree-level ef- 
fective potential: due to C and CP conservation it is a block-diagonal 6x6 
matrix containing the 3x3 neutral scalars matrix, the 2x2 pseudoscalar ma- 
trix and the charged scalar mass. The elements of these matrices are given 



Ay^^-r^(0,T)~;^ 



F,B (27r) 2 




m ~ 



ITeV 



by: 



a)neutral scalars: 
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in the basis {ReHi, ReH2, ReN) we have 



^lii = 2ml^ + 2\\W+\Hl\^) 



Mj22 = 2m'+2X\\N\' + \Hl\') 



MI33 = 2ml + 2\''{\Hl\'' + \Hl\'' + \H+\^ + \H-\^)+Ak''[2>{ReNf + {ImNf] 
- AkAkReN - A\kRe{H^H2); 



Ml 



s 12 



= [4A2 - 2(c/2 + g'^)]ReHfReHl - AXA^ReN 

- AXk [(ReNf - (ImNf] - 4 (A^ - ^) Re{H+H-); 



Ms 13 = AX^ReH^ReN - AXkRe{H^N*) - 2XAxReH^; 



Ar*^ 



ro. 



Ms 23 = AX'ReH^ReN - AXkRe{H^N*) - 2XAxReH'l; 



b)pseudoscalars: 

in the basis (^40, ImN) 



p 11 



sin^ P 



;(i4) 
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Mp22 = m% + +2\WHi\^ + \H2\^) 

+ 4A;2[2|iV|2 - Re{N^)] + AXkRe{HiH2) + AkAkReN; (15) 

Mp 12 = sin (3 [AX^ImNImHl - 4XkRe{N*H^) + 2XAxReH^] 

+ cos p{'iX'^ImNImH^ - 4XkRe{N*H^) + 2XAxReH^} ; 



c) charged scalar: 



ml 



cos^ P 



mjj^ + X^{\N\^ + \H-\^) + 



£^-^1 0-0,2 

A l^ll 



+ 



9' + 9"nrr0^2 



{\H!i\' + 2\H^ 



+ sin (5 cos (5 



+ ^{H^H^) + XkN^ + XAxN* + h.c. 



+ sin^/S 



+ 



^2 + g'^ 



{\El\^ + 2\E-\'-\E^\') 



Next we come to the SUSY particles. Among the squarks and leptons 
we should consider only the stop-quarks t/^^ij, the others having negligible 
Yukawa couplings. However fj, r masses (9) are dominated by ttiq, which at 
the electroweak scale is much greater than M^y (its renormalization group 
equation is dominated by the strong coupling [19]) and then their contribu- 
tion is Boltzmann-suppressed. 

Charginos have a 2x2 mass matrix (in the basis (-f^~, W^~)) 



2 _ / AiV g^E 



(16) 
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where M2 is the gaugino direct mass term. Assuming 

M2 f5i 771 » Mw ~ Ax 

we get 

which shows that the hghtest chargino may be in thermodynamic equihbrium 
and contributes to (13). The LEP bound on the lowest chargino mass [23] 

m^± > AbGeV 

will be imposed to constrain the parameter space for A and x. The 5x5 
neutralinos mass matrix, Ml, has two heavy (of order Mi^2) Boltzmann sup- 
pressed eigenvalues. With the choice of the parameters that will be illustrated 
in the following, the off-diagonal terms Mn 12^, 13^ and 23^, are negli- 
gible in comparison with the three diagonal ones, n^. Mat 22^ and 33^, 
so that the approximate masses for the three lightest neutralinos are given 
by the following expressions 

m^,' ^ \\\N\' + \H',\') + {g' + g")\Hl\'- 
m^,' ^ \\\N\' + \Hl\^) + {g^ + g")\HT. 
mm' ~ X\\H',\' + \H^,\')+4e\N\\ 

Now we are ready to write the one-loop, T ^ 0, correction in the high 
temperature limit. It is given by 

AV^^o°P{(p,T) = — [TrM| + TrM|, + 2m^ + 6M^ + 3M| 
+ Gm'^t + 2mlhi + + m%2 + m%s]' , 

where [■■■]' means that the sum runs on the light masses only. As the field- 
dependent mass of a given particle becomes much greater than the tempera- 
ture, the contribution of that particle to the previous expression is switched 
off. Numerical evaluations [9] of the integrals in (11, 12) show that we can 
approximate them at better than ten percent using expression (13) for masses 
up to m/T ~ 2. 
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4. Determination of the critical temperature. 



We defined the critical temperature Tc as that value of T at which the 
origin of the field space becomes a saddle point for the effective potential. 
In fact, it is well known [9,24] that the transition proceeds by quantum tun- 
nelling and is completed at a temperature higher than Tc, however we are 
now interested in determining the order of the transition, which is parame- 
terized by the vacuum expectation value at Tc {v{Tc) = means a second 
order phase transition, v{Tc) 7^ a first order one). 

It is easy to check that at high temperature symmetry is restored, in the 
sense that the origin is a minimum of the effective potential. 
The critical temperature is then defined by the condition 



where the effective mass matrix Mj is given by the second derivatives of the 
full one-loop potential at finite temperature with respect to the scalar fields. 
In the origin of the field space the effective mass matrix is approximatively 
diagonal in the basis {ReHi, ReH2, ReN) 



det \Mf{Tc) 



J 0i=O 



= 0, 



(17) 



2 










\ 



where 




9A2 + 12A;2 + 3AA;sin2/3, 



and 



A, 




~ 2 
m 



rad 



87r2 



,2 -I- mf 
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Neglecting the off-diagonal term, the critical temperature is given by the 
highest among T^^, T^^ and where 

Tf ^ JJnilAHl^ (i^H,,H,,N). (19) 

Obviously the effective potential can become ffat only along those directions 
corresponding to negative soft masses. Negative values of may however 
be dangerous for CP conservation, which we have assumed as a 'boundary 
condition' for the vacuum at T=0; if we look at the effective mass matrix for 
pseudoscalars in the origin of field space, and in the basis (^4°, ImN), we get 



M 

where 



m 



%^ sin^ /? + mjj^ cos2 (3 + ^Pao 



/ io 

= + I2k^ - 3AA; sin 2/?. 



Then, at the temperature 

p2 _ 2Am% s2 

-TAT 

the potential becomes flat, in the origin, along the direction ImN, while 
it remains convex in the direction ReN. So, if at that temperature the 
transition in the scalar fields has not yet taken place then the imaginary part 
of the fields acquires a vacuum expectation value, and, when the temperature 
goes to zero, the vacuum is no more (pi = (px = (see eq.(5)). 

Due to the heavy top [19], can run to negative values at low energy, 
whereas, if A is not too large, m|^^ remains positive. It seems also plausible 
that, for small A and k also mj^ remains positive, although a detailed reso- 
lution of the RGB's for the soft parameters would be necessary in order to 
clarify this point. 

^Note [20] that in {ReHi, ReH2, ReN) = ivi,V2,x) the pomt {A° , ImN) = (0,0) is 
always a minimum, as the pseudoscalar mass matrix has only positive eigenvalues. 
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However, in what follows, we will constrain the parameter space by means 
of eq. (6), requiring that the only negative soft mass is m^^- 
The critical temperature will then be given, in any case, by 



where the is due to the fact that we have neglected the off-diagonal term 
in (16). Anyway, numerical results will be obtained by solving the condition 
(17) with the complete matrix. 

5. Minimization of the potential 



Once we have ensured the correct pattern of symmetry breaking, i.e. no 
VEV's for ImN, and C+ (this last condition is guaranteed by eq. (7)) we 
can study the minimization of the potential, at the critical temperature, in 
the scalar directions only. 

Before illustrating the numerical results of this minimization let us give 
an intuitive view of what actually happens, namely, of the crucial role played 
by the trilinear soft SUSY breaking parameters A\ and Ak in the formation 
of a minimum of the effective potential for large values of the fields (~ Tc). 
We define the polar coordinates 

ReH^ = Y cos a cos 
< i?ei?2 = i^cosQ! sin/3^; 
ReN = Ysina. 

The full effective potential may then be written as the sum of terms of dif- 
ferent order in Y = Y/T 

V=^ = aY + bY^ + cY^ + dY^ + ---, (20) 

where the ellipses indicate terms coming from AV^~q°°^, which, in our ap- 
proximation, are different from zero in the direction if^ only and, due to the 



re « Tl 



\ 
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renormalization condition (10), will be small as long as we search a minimum 
at Tc with if 2 t^e same order of f 2 . 
The coefficients in (20) are given by: 

i) linear terms: 



a = - sm 2p sm a-^; 



ii) quadratic terms: 
1 



6 = 



2^2 



cos^ a{niH^ cos^ + mjj^ sin^ 0^') + mj^ sin^ « + ^(- • •) 



2 /3T\ 



2i2 



24 



here we do not write explicitly the terms coming from finite temperature cor- 
rections; for the present purposes what matters is that the complete quadratic 
term is positive or zero at T = Tc in any direction; 

iii) cubic terms: 

1 / 2 \ 

c— — — sin q; iXAx cos^ a sin 2/3^ + -kAk sin^ aj ; 

iv) quartic terms: 

d = cos^ a ^sin^ + ^ ^^^^ ^^"^ 2/3^^ + /c^ sin^ a 



■sin^2Q;sin2/3^ + 



cos'* CKCos^ 2/3^. 



At large values of the fields, > a/c{^ (that is F > Tq) the potential 

has a minimum at 

~ 3c 

Y PS ; 

4d' 

the condition for this to be a global minimum is (neglecting the quadratic 
term) 

-3 sin^ a iXA\ cos^ a sin 23'^ + ^kAu sin^ a] 

^ . ^> 1 



ad? 



Axsm2(3T^\^ ^ ^^^^ 

Assuming that eq.(21) is satisfied, the question is now to determine the 
direction of the vacuum in the field space. Since the only negative term 
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appearing in (20) is the cubic one, and it is proportional to sin a , there will 
be always a i?eA^-component in the direction of the vacuum. 

However, for large Ax the l.h.s. of (21) decreases for sin a = 1 (it is 
the effect of the linear term that raises the potential in that direction), so 
that the sin a directions become favorite, and an electroweak first-order 
(nonzero VEV's for ReHi, ReH2) phase transition takes place. 

6. Numerical results. 

Contrary to the case of the Standard Model, the constraint (2) on the 
vacuum expectation value of the fields at the critical temperature cannot be 
written as a simple bound on the masses of the scalars in this model. In 
fact, approximate analytic expressions for these masses can be obtained [20] 
in the limiting cases x Vi, ^2 or a; <C fi,f2, but we are here interested 
in values of x of the same order of the electroweak VEV's Vi,V2 0. So we 
are induced to perform a numerical investigation on the parameter space in 
order to find the region in which the bound (2) is satisfied, and then look at 
the corresponding values for the mass of the lightest scalar. 
In addition to the six tree-level parameters, we have to consider also those 
appearing in the one-loop corrections, namely, the top Yukawa coupling ht, 
the stop soft SUSY breaking mass rhq, and the gauginos direct masses Mi 
and M2 (14). We will assume a common value for rhn, Mf and M| 



Then we constrain ht, X, and k by means of the renormalization group anal- 
ysis performed by the authors of ref [25]. Requiring that the coupling A 
remains perturbative up to a large (say 10^^ GeV) scale, they find, at Mz 



^One of the original reasons for the introduction of this model was the so called /i- 
problem, whose solution requires Xx ^ 0{Mw) [20]. 



ml = Ml = Ml = m^ = ITeV. 




2Ml 



(22) 



and 
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ht{Mz 



0.97 



(23) 



Moreover, RGB's for the couplings have the fixed ratio point (neglecting g, g' 
with respect to ht) 

We will fix A at its upper value as given by (22) and ht and k according to 
(23), (24). 

For fixed M2 and A the experimental limit [23] on the chargino mass gives 
a lower bound on x, approximatively 



A^GeV 
A~' 



(25) 



We then fix tan^ to some typical value (we will use tan(5 — 2, 10) and require 
(see the discussion on the critical temperature) 



> 0, w?H, < 0. 



This, in turn, implies 
where 



(26) 



A mm 



A max 



1 

X 
1 
X 



Xcotg(3 ( + ^ 
XtgP (x^ + ^ ) - 



+ — sin 23 - kx'^ 
4 



Xv 



sin 2/3 — kx"^ 



or, equivalently 



A 



sin^/3 



X + — \ <mi < 



cos^ (3 \ 4 



where is the mass of the charged scalar. Eq. (26) automatically ensures 
that eq. (8) is satisfied, i.e., that the charged fields do not acquire any VEV. 
Finally, the requirement m% > implies a lower bound on Aj. 



Ak > 2kx + 



Xv^ 
kx 



X-sm2p(^ + k 



(27) 
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For any triplet {x,Ax,Ak) satisfying (25), (26), and (27), we minimize 
numerically the effective potential at zero temperature and verify that the 
vacuum is at {ReH^, i?ei?2 , ^eA^) = i, V2, x). We then calculate the critical 
temperature Tq according to the definition (17), and, finally, minimize the 
effective potential at Tq- 

In Fig. 1 wc report v{Tc)/Tc against the lightest scalar (^/isi) mass for all 
the points that passed the above mentioned selections. As we see there are a 
lot of points corresponding to v{Tc)/Tc > 1.3 and mgj > &QGeV. Note that 
the maximum allowed values for m^^ reach 170^61^. 

In Figs. 2a, 2b we plot the allowed region (i.e. the points corresponding 
to v{Tc)/Tc > 1.3) in the ms^ — rric plane, where rric is the charged scalar 
mass; x is fixed at 300 GeV, while the upper and lower limits for nic are 
given by (26). 

Finally, in Fig. 3a, 3b we recover the result we got from the naive discus- 
sion on the minimization of the potential, namely, the crucial role played by 
A\ in the determination of the direction of the transition. The dashed lines 
represent the values of v{Tc), while the continuos ones those of Tc- 

In conclusion, we have showed explicitly that minimal extensions of the 
MSSM allow to easily circumvent the Higgs mass problem arising from the 
study of the electroweak phase transition in the Standard Model and in the 
MSSM. The main role is played by the trilinear soft SUSY breaking terms, 
which are present in any supersymmetric model with extra SU(2)-singlet su- 
perfields. Radiative effects are also important as, due to the heavy top, they 
force the direction of the transition (i.e. the fiat direction of the potential at 
Tc) to be always ReH2. 
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Figure Captions. 



Fig. 1 v{Tc)/Tc ratio vs. the lightest scalar mass, for all the allowed region 
for the parameters A, k, x, ^4^, and Ak (see text), and tgj3 — 2. 



Fig. 2a Upper and lower bounds for the lightest scalar mass coming from 
the various conditions imposed on parameter space (see text) plus the re- 
quirement that v{Tc)/Tc) > 1/3. x has been fixed to 300 GeV and tgP — 2. 
Fig. 2b The same as in fig 2a with tgP — 10. 



Fig. 3a v{Tc) and Tc dependence from Ax, for x = SOOGeV, Ak = 190GeV 

and tgf] = 2. 

Fig. 3b Same as in Fig. 3a for x = AQQGeV, Ak = SOQGeV and tg(3 = 10. 
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